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We investigate some aspects of the self-consistency in the Dyson-Schwinger approach to both the 
QED and the self-interacting scalar field theories. We prove that the set of the Dyson-Schwinger 
equations, together with the Green-Ward-Takahashi identity, is equivalent to the analogous set of 
integral equations studied in condensed matter, namely many-body perturbation theory, where it is 
solved self-consistently and iteratively. In this framework, we compute the non-perturbative solution 
of the gap equation for the self-interacting scalar field theory. 
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I. INTRODUCTION 



Nonperturbative approaches to quantum field theory (QFT) allow for a better understanding of some general prop- 
erties of the exact scattering amplitudes. They historically pursued after the formulation of quantum electrodynamics 
(QED) to analyse the asymptotic behaviour of renormalized operators at short distance (see, for instance, Ref. [l|). 

Beyond that, they give considerable informations in detailed studies of the structure of higher order approximations 
in theories where the perturbation expansion, along with a strong coupling, fails in analysing the short distance 
behaviour of the relevant theory operators. This approach (referred to as Dyson-Schwinger) has been historically 
introduced by a number of authors 

Since then, it has inspired plenty of papers on the subject (see, for instance, Refs. 043)- Among them we quote 



those dealing with non-abelian theories jl0l - ll5l |. such as quantum chromodynamics (QCD), where both the aymptotic 
freedom (large coupling values at low energy) and the confinement of the Lagrangian fields within the asymptotic 
ones intrinsically require a phcnomcnological analysis beyond the naive perturbative expansion. In these processes 
the short-distance contributions can be computed, to some extent, by using the factorization approximation for the 
local operators in the effective Hamiltonian, as argued by Bjorken on the basis of color-transparency (l6j | . 

Large coupling and field confinement also characterize solid state physics phenomenology, where a screened many- 
body interaction occurs at about the Fermi energy among a large number of electrons embedded in the crystal lattice. 
In condensed matter singularities appear in the perturbation series of the correlation energy of a fully degenerate 
Fermi-Dirac system with Coulomb interaction, otherwise named homogeneous electron gas (HEG): in Ref. [17| the 
large divergent logarithms are resummed according to the general scheme introduced by Feynman in QFT ja]. In 
spite of its simplicity, HEG is fundamental for computing the correlation energy of a wide class of complex systems 
within the so called local-density approximation [18j |. 




Also within this context, a Green's function theory approach, called (improperly) many-body perturbation theory 
(MBPT), based on a formalism of second quantization of operators, has been considered [l9|, [2y]. The fundamental 
degree of freedom is the Green's function or propagator, which represents the probability amplitude for the propagation 
of an electron. As in any other QFT, the many-body system can be expanded in perturbation theory, with the coupling 
being the many-body interaction term. The Green's function (as well as any other quantity of the theory, such as the 
self-energy or the polarization) can be calculated at a given order of perturbation theory. A Feynman diagrammatic 
analysis is hence possible. The theory at the first order is equivalent to the Hartree-Fock theory. However the coupling 
is not small (compared, for example, to the electron-ion interaction) and the expansion does not converge. The second 
order is not necessarily smaller than the first. Hence one needs to resort to more complicated methods to solve the 
theory. Beyond partial resummations of diagrams to all orders, iterative methods have been preferred. 

Historically a Dyson-Schwinger approach has been introduced to account for such a non-perturbative phenomenol- 
ogy and to compute optical and electronic properties of complex systems by means of iterative schemes. The alternative 
formulation in terms of functional derivatives reduces the many-body problem to the solution of a coupled set of non- 
linear integral equations, whose characteristic feature is that, besides the single-particle Green's function, a whole 
hierarchy of equations involving higher order Green's functions is generated. A truncation of Feynman diagrammatics 
corresponds to the custom of replacing this hierarchy of equations by another coupled set of nonlinear equations 
connecting the single-particle Green's function to the mass, polarization and vertex operators, often referred to as 
the Hedin's equations in condensed matter [2l|, to be solved iteratively until self-consistency is achieved^]. Here the 
higher order Green's functions dependence of the relevant quantities is recasted within the functional derivative of the 
mass operator with respect to the fermion single-particle Green's function. The latter quantity is proved to equate 
the Bethe-Salpeter kernel for the two-body into two-body rescattering. Both the diagrammatic truncation and the 
nonlinear coupled set of equations correspond to a nonperturbative approximation to the solution of the many-body 
problem. 

So far, nobody has solved the Hedin's equations for a real system, since the nonlinearity of the equation involving the 
Bethe-Salpeter kernel is computationally demanding. Approximations are required to simplify the problem. Among 
the most widely used computation schemes it is worth mentioning the so called GW approximation, where the vertex 
operator is simplified in the self-energy evaluation at the beginning, and the Bethe-Salpeter equation accounts for 
the vertex corrections within this approach (see, for instance, the reviews in [IH |22|). 

The analogies between these two fields, spanning ultrahigh and ultralow energies, have been cross fertilizing method- 
ologies and approaches. Along with that, this paper aims at spotting: i) the equivalence between the Hedin's equations 
set and the Dyson-Schwinger approach to QFT, where the role played by the underlying gauge symmetry is crucial 



1 The approach described here is akin to the truncation procedure of the infinite tower of Dyson-Schwinger equations in particle physics 
(see, for instance, Ref. 0). 

2 The name G W stands for the product of the G Green's function and the W dressed Coulomb interaction while computing the mass 
operator YGW , being T the vertex operator as it appears in the Lagrangian (T w I when magnetic and relativistic effects are neglected). 
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in relating the Bethc-Salpctcr kernel to the functional derivative of the mass operator with respect to the fermion 
Green's function (throughout the Green- Ward-Takahashi identity [23T - [25l | ) ; ii) the iterative scheme solution of the 
Dyson-Schwinger equations, in the spirit of the condensed matter methods, by means of a S-matrix unitarity in- 
spired ansatz on the mass operator, accounting for the nonlinearity of the present approach. Both points should 
be regarded as aspects of the self-consistency of the Dyson-Schwinger approach to the many-body quantum field 
theory, with regard to the actual theory investigated, in the context of the aforementioned S-matrix unitarity inspired 
picture. Indeed, in QFT the S-matrix satisfies the unitarity condition, actually more fundamental than the concept 
of Hamiltonian and wave functions (see, for instance, Refs. [U I26l - l28j and references therein). 

In that respect, we point out two examples (QED and the self-interacting scalar field theory), different on the 
physical basis but akin as to the solution scheme. The QED example is intended to formally bridge the condensed 
matter scenario to its underlying fundamental theory: throughout the paper, a QFT analysis to the topics has been 
preferred to the functional one, in order to spot this idea. A functional approach is, of course, viable. On the other 
hand, the gap equation arising in the self-interacting scalar field theory is probably the simplest example to show how 
to implement the iterative scheme to the numerical solution of the Dyson-Schwinger equations. 

The plan of the paper is as follows. In the next section we discuss the general formalism of Dyson-Schwinger 
approach in QED. In section 3, the Green- Ward-Takahashi identity is proved to be equivalent to the definition of the 
Bcthe-Salpeter kernel as the functional derivative of the mass operator with respect to the Green's function. Section 
4 is devoted to the calculation of the Dyson-Schwinger equations for the self-interacting scalar field theory. Finally 
in section 5 we compute the non-perturbative solution to the gap equation arising for the latter case and draw our 
conclusions. 



II. THE SET OF DYSON-SCHWINGER EQUATIONS IN QED 



It is worth to recall the Dyson-Schwinger equations for QED. Hereafter, latin letters shall refer to noninteracting 
quantities, while calligraphic symbols represent interacting ones. Thus, ip represents the free fermionic field, while ^ 
is the exact fermionic field; analogously for the photonic field A (A). Accordingly, G (Q) is the free (exact) fermionic 
Green's function, while D (V) corresponds to the free (exact) photon propagator. The electromagnetic current j is 

defined according to Ref. HH, i.e. ^{x) d = ^>(x)^^>(x). They read: 

G%{x,y) de d- -i{O\T(r{x)^ {y))\O) , (i 9 * - m I) #n) = , (1) 

D»»(x,v) +i(Q\T (A»(x)A^(y)) |0) , {VP - d a ff*gT)A v (x) = , (2) 

G<${x,y) de J- -i{Q\T(* a (x)^ p {y))\Q) , (i X - m X) *(*) = eJL{x)*{x) , (3) 

V^(x,y) d ^ +i(0\T(A ft (x)A^(y))\0) , (0"0* - a W) A v {x) = ~4ne f(x) , (4) 

(0\T{A»(x)* a (y)*p(z)) |0) de d- e JJJ d 4 x'd 4 y'd 4 z' \z,z') ry'(z\y\x') G$tf ,y) V»,(x',x)} , (5) 

M$(x,y) d ^ G/ a (x,y)-gp la (x,y) , V^{x,y) = 7 ' D^^x, y) - V~ l ^{x,y) , (6) 

where a sum over repeated indices is understood. The relevant equation for M. operator is obtained by applying 
= ' 8 4 (x — y) (i $ y — mX)g J to Q operator defined in Eq. ((3]) and by using Eq. ([5|) for the vertex Tq 

It reads: 

M1{x,y) = -ie 2 ff d 4 x'd 4 y' 7 ^ G$(x,x>) T^x' , y , y') V^{y',x) ; (7) 

analogously for V operator one gets: 

V^(x,y) =47rze 2 JJ d 4 x'd 4 y' 7 J ^' (*,*') T^{x',y',y) g*{y',x) . (8) 

Beside these, another equation is needed for the vertex F: 
iTp^ix, y, z) d = i^5 4 {x - z)6 4 (y -z) + 



3 Details on this calculation are carefully reported in Ref. [29ll , 
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/ dVdVdV'A" £#W) r«>(x',y',z) g^{ V ',y") K^(y,y",x,x") , (9) 

where K, the Bethe-Salpeter kernel, accounts for all the possible contributions coming from the fermion-fermion 
rescattering, except for those already embodied within the exact propagators G,T> and vertex 10. 
Equations (|6I7I8I9[) are also reported in a diagrammatic representation: 







M 






M 


def. 
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def. 



def. 



• + 




(3) 



As a matter of fact, the set of five integral equations depicted above depends on six unknowns, i.e. Q,T>,V,A4,T 
and K. Thus, an additional equation is needed to close the system and to approach a solution. Two schemes have 
been assessed to accomplish with such an issue. In condensed matter the Bethe-Salpeter kernel is written as the 
functional derivative of the mass operator with respect to the Green's function: K = 5M./8Q [lj| r2ll. I22I . In QFT, 
the same definition formally holds. Indeed, in the framework of Cornwall-Jackiw-Tomboulis approach |30j|. it is known 
that the second functional derivative of the action with respect to the (bilocal) field equals the Bethe-Salpeter kernel 
[3lj . On the other hand, in Ref. 11 1 such a second functional derivative is proved to be the functional derivative 
of the mass operator with respect to the fermion Green's function. While such a relation is crucial for the self- 
consistency of the coupled set of nonlinear equations, in order to avoid a hierarchy of higher order Green's functions, 
a functional derivative is computationally deman ding and a diagrammatic, although non-perturbative, expansion of 
the Bethe-Salpeter kernel is usually preferred!! [7l4l0l Ii~2r4ll| . 



III. THE GREEN-WARD-TAKAHASHI IDENTITY. 

Under a small phase (gauge) shift on the ^ operators, Sx, the G{x, x') operator is shifted by ieQ(x, x) [S\(x) — 5x{x')), 
according to the definition of Eq. On the other hand, 8Q(x,x') can be directly computed as the coupling of the 
(small) gauge field —d^Sx to the fermionic current By equating the former and the latter quantities, one gets the 
Green's equation [24j either in its integral formulation, i.e. 



ieg$(x,x')[6 X (x) - 5 X (x')] 



d 4 yd 4 y'd A z g%(x,y) d^Tf^y , y' , z) \y\x') S X (z 



(10) 



or in its differential one — ieQ 1 (x 1 x') [8x{x) — 8x{x')\ = e J d 4 z d^T^x, x' , z) Sx(z). We will need both forms in the 

sequel. The Ward-Takahashi identity corresponds to the soft photon limit in momentum space (23l . [25j . 

To accomplish with the task of completing the set of Dyson-Schwinger equations, we replace Q~ l in the differential 
form of Green- Ward-Takahashi equation, according to the mass operator definition of Eq. (|6]), and we notice that 
the non- interacting propagator G is gauge invariant, therefore 



iG- 1 (x,x')[6 X (x)-6 X (x')} 



d 4 z 8 X (z) d» [l^ 4 {x - z)8 4 (x' ~ z)} 







(11) 



4 In a functional approach they are named irreducible diagrams. 

5 In particle physics, this approach is usually referred to as truncation of Dyson-Schwinger equations. 

6 Polarization does not affect the gauge field, since the former is a transverse tensor while the latter is longitudinal. Once again, details 
can be found in Ref. 1231. 
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Thus: 

ieM(x,x')[6 X (x)-6 X (x')] = e J d 4 z 6 X (z) d% [-^S 4 (x - z)5 A {x' - z) + T ll (x,x',z)] = 
= e j d 4 z S X {z) JJJJ d 4 x'd 4 y'd 4 x"d 4 x"Q{x',x") 0JT„(x", y", z) G(y" ,y') K(y', x', x', x) = 

d 4 x'd 4 y' Sg(x\y')K(y\x',x',x) , (12) 



where in the latter equations we used the vertex definition of Eq. © and the integral form of Green- Ward-Takahashi 
relation of Eq. (fTDj) . It can be easily seen that the l.h.s. of Eq. (|T2l is the shift of the mass operator under a gauge 
transforrrQ, i.e. SA4(x, x'). 
In conclusion 

SM(x,x') = JJ d 4 x'd 4 y' 6G(x',y') K(y' , x' , x' , x) . (13) 

The meaning of this equation is twofold. It states that: i) the gauge shift on the mass operator is linearly dependent on 
the gauge shift of the fermion propagator throughout the Bethe-Salpeter kernel; ii) for the same reason, the Bethe- 
Salpetcr kernel is gauge invariant. Moreover, the equation written above completes the set of Dyson-Schwinger 
equations and makes finding a solution a viable problem. 

IV. THE SET OF DYSON-SCHWINGER EQUATIONS IN THE A(0*0) 2 THEORY. 

The set of Dyson-Schwinger equations for self-interacting scalar field theory is simpler than the analogous in QED. 
Hereafter <p (0) represents the free (exact) bosonic field and, accordingly, D (T>) is the free (exact) boson Green's 
function. It reads: 

D(x,y) d ^ +i(0\T(cp(x)<p*(y))\0) , (□, + m 2 ) ip{x) - , (14) 

V(x,y) d ^ +i(0\T(<f ) (x)cj ) *(y))\0) , K+m 2 )^)^^*^ 2 ^) , (15) 

(O|T(«^(y)0*(z)0*M) |0> = 7 ' -iX x 

d 4 x'd 4 y'd 4 z'd 4 w' {V{w, w') V{z, z') IL(w' , z', y', x') V(y', y) V(x' , x)} (16) 

M(x,y) d ^ D- 1 (x,y)-V- 1 (x,y) , (17) 

d d 

being O x =' — — — — . The latter and the former sets are akin as to the derivation. Let us consider, for instance, the 

OX^ OXfj, 

mass operator equation. By applying D~ x to T> (eqs. I14I15|) one gets: 

.A 
'2\ 



d 4 y D-\x,y) V{y,x') = S 4 (x - x') + i± (0|T (<j>{x) 2 <t>* {x)<t>* {x')) |0) . (18) 



Finally, by using Eq. (|16p . the mass operator is computed: 

M(x,x>) = ^ JJJ ' d 4 y'd 4 y"d 4 x" V{x,f)V{x,x'')U{y'',x'', y ',x')V{y',x) . (19) 

The vertex equation introduces a Bethe-Salpeter kernel for the two-body rescattering: 
n(y', x', y, x) = 5 4 (y' - x')S 4 (y' - y)S 4 (y' - x) + 

+ {HI d 4 x'd 4 x"d 4 y'd 4 y" B(y', x' , y", x") V{y" , y') V{x" , x') K(y', x', y, x) . (20) 



7 Indeed under a gauge shift: 5 ^ J d 4 y Q(x,y)Q 1 (y,x') j = 0. By expanding the latter gauge variation and by noticing that the 
non-interacting propagator G (unlike Q and Ai) is gauge invariant, we achieve the result. 
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Although the self-interacting scalar field theory is not a gauge theory, the gauge shifts of mass operator and 
propagator (under a phase shift of <p operator) are related in the same fashion of QED cas^J 



(21) 



The latter equation completes the set of Dyson-Schwinger ones , Eqs. (fT71 fT9l l20|) for the self-interacting scalar field 
theory. They can be recapitulated in the following diagrammatics where S [• • •] means the gauge shift of the quantity 
within brackets: 



© 




def. 



• + 




© 



lef. 



5[@\ 



K 



(1) 



V. RESULTS AND DISCUSSION 

Hereafter we shall assume that the space-time is homogeneous. Therefore the energy-momentum conservation law 
holds in momentum space, where the relevant quantities defined in the previous section simplify accordingly. Thus 
the set of Dyson-Schwinger equations reduces to: 

M{p) de d- D-\p)-T>-\p) , 



M(P) 



n(p + cr,cj - <r,p, q) 



def. 



2! .// (2tt) 4 (2tt) 
d 4 k 



— £>(c) T>(rf) n(-cr, -ri,p - a - rj, -p) V(a + t]-p) , 



(2tt) 4 



U(p + a,q - <7,p + k,q - k) V(q - k) V(p + k) K(p + k,q — k,p, q) , 



M(p + q) -M(p) 



d 4 k 



V{k)-V{k~q) K(k-q,p + q,k,p) 



(22) 



For the sake of clarity, we prove the last equation in (f22|) . We start from the last equation in (fT2"|) for the self- interacting 
scalar field theory: 



SM(x,x') 
SM{k,k') 



d 4 x'd 4 y' 8V(x',y') K(y',x',x',x) 
d 4 k d 4 k' 



(momentum space) 
6V(k, k') K(k', k', k, k) <=> (5% is the gauge shift) <=^> 

/ d 4 q/(2-K) 4 6x(q) X 



(2tt) 4 (2tt) 4 

d 4 q/(2n) 4 Sx(q) \M(k + q, k') - M{k, k' - q) 
d 4 k d 4 k' 



T)(k + q, k') — T>(k, k' — q) K(fc', k', k, k) (homogeneous space-time^) 

d 4 k 



(2tt) 4 (2tt) 4 

i e 6x(k' - k) \M(k') - M{k) = i e S x (k - k') 
x K(fc + k - k', k', k, k) 



(2n) 4 



V{k + k-k')-V{k) 



(23) 



8 For a strictly neutral </> field, i. e. <j>* = (f>, an additional (Schwinger) term arises. The absence of a global gauge symmetry for this case 
has been also pointed out in [32l |. 

9 For instance, the propagator in the momentum space reads V(k,k'); it reduces to (2n) 4 t> (k) 8 4 {k - k'). On the same footing, the 
gauge shift ST>(p,p') = ie J d 4 q/(2w) 4 &x{q) \j}{p + q,p') — T>(p,p' — g)j is shortened to &x(p' ~ p) \^(p') ~ ■ Analogously for 
the mass operator. The Bethe-Salpeter kernel is also simplified: K.(p,q,p' ,q') = (2tt) 4 S 4 (p + q — p' — q') K(p, q, p' , q'). 
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after a variable shift (k' —> p + q, k — > p, k —> k), and accounting for the ^Y-oddnesJ^l. the proof is completed. 

The S- matrix unitarity naturally leads to physical conditions for the propagator structure, i.e. a pole term times 
a multiplicative constant (renormalization residue) [2j|. Within this context, although not necdecO, we shall in- 
troduce the following ansatz on the mass operator: A4(p ) = — aD~ 1 (p 2 ) + /i 2 , being a a (real) constant to be 
determined, while fi is a mass scale (renormalization point). By this choice, the interacting propagator reads 



V~ x {p 2 ) = (1 + a) D-^p 2 



1 



M( P 2 ) 



and the pole is located at D 



\Ppole 



/'" 



i 



More- 



1 + a 

over, by using this S-matrix unitarity inspired ansatz, the Bethc-Salpctcr kernel K is promptly evaluated by means 
of the last equation in (f2"2")l : 

a 

q,p- 



K{k 



q,k,p) 



1 



V^ik-q) V^ik) (2n) 4 5 4 (p- k + q) 



(24) 



by which the vertex II is shortened to (1 
of (EH). 

Finally, the mass operator equation in 



a) I (vertex renormalization) by direct substitution in the third equation 



M{ P 2 ) = (1 + a) 



accounts for the self-consistency of the Dyson-Schwinger approach 
A 2 



2! 



d A x exp (i px) \D(x)\ 



(25) 



where T>(x) is the interacting propagator in space-time (configuration space), while the implicit dependence upon a in 
both T> is understood* 2 !. Eq. (|25| describes a setting sun Feynman diagram: it provides either the broken symmetry 
solution or the symmetric one. To accomplish with the task, the mass operator of Eq. (|25[) satisfies a dispersion relation 

( +oo dt Im (M(i) N 
(DR) with one subtraction |33H35| (in p variable Ji3 Msub{p 2 ) = p 2 ' 



version of Eq. (1251) on the mass pole, resembling the mass/gap equation of [36l 13 



P 2 ) 



We evaluate the subtracted 



A 2 * 



2!(2tt) 4 



+oo 



dx 



J\ (mi) 



In euclidcan space it reads [38| : 
K^mxf x J-(x,A) (26) 



where J\ is the first order Bessel function of first kind, Xr = A/(l + a) and J-{x, A) is an ultraviolet cut-off factor |36|. 
The integral of Eq. (|26|) has to be regularized at some point (1/A) since it exhibits a (logaritmic) divergence for small 
x (A — > oo) (the variable x behaves like the inverse four- momentum p, thus at the origin x ~ 1/A). 

After rescaling the variable x — > £/m and normalizing the parameter m to the cut-off A (m — > rj A), Eq. (|26p 
reduces tcF^I: 



2 2 
V =V 



2!(2tt) 4 



4-oo 



(27) 



For vanishing 77 the mild divergence of the integral on the r.h.s. of Eq. (f2"T|) is compensated by the n 2 term and the 
trivial solution (77 = 0) is achieved. 

The importance of Eq. (|27p is twofold. On one hand, a mapping between \r and r\ can be obtained beyond the 
perturbative regime of the coupling constant, helpful while investigating the cohexistence of the triviality (namely 
Xjj — > when A — > 00) and the spontaneous symmetry breaking (SSB) of the theory [39| (crucial, for instance, for 
the self-consistency of the Standard Model in particle physics |40jV 

On the other hand, a renormalization group equation (RGE) for the 77-evolution of the coupling constant \r can 
be promptly assessed from Eq. ([27]) : 



A|(*7) 



Ajjfao) 



77—^0, f]o/r] fixed 



A|(r?o) 



2lf27r) 4 Z?) 7 ^ X ^C 3 -) 



i+ 2!(4tt)1 ln VTj 



(28) 



10 The phase (gauge) shift oddness stems from the Green's function definition : T>(x',x) = — T>(x,x')* . 

11 In principle, the set of equations accounts for an even number of unknowns (M, T>, Tl and K): it can be numerically solved 
sclf-consistently. 

12 In euclidcan four-dimensional space the interacting propagator, together with the mass operator ansatz, reads T>(x, m) = — — , 

1 + a (27r) 2 x 2 

m = , x = \ /x^x^, being K\ the first order modified Bessel function of second kind, while M(p ) = cup + fi . 

y/l + a v 

- -9 f + °° dt Im (A/(t)) 

13 Otherwise M su b{p ) = M(jr) — M(0), being M{p ) = / 5-^- the mass unsubtracted DR. 

J Q tt t-p 2 

Hereafter we implement ^(x, A) = 9 [x — 1/A], being 6 the Hcavisidc function. 
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being f2(x) the integrand function of Eq. ([27]) . While the latter equation corresponds to the perturbative region for 
the coupling A^, the former one holds in non-perturbative regime too. 

Beyond the trivial solution (ry = 0), corresponding to the unbroken symmetry [4l|, Eq. ([27]) provides with the broken 
symmetry solution, once the renormalized coupling \r is replaced according to the SSB mechanism for the theory 

investigated. Indeed both m 2 and Xji depend upon the renormalization constant a in the same fashion [ex 

V 1 + a 

2 f_l 2 

thus, by eliminating it, we get the formula m = — A/?, accounting for the not vanishing vacuum expectation value of 

A 



the field which minimizes the Hamiltonian ^(O|0*0|O) = v ^ 0, — = — ) . It reads a typical fixed point problem 

for the parameter rh ( rh d =' — , A d = ' — | [42l, III : 

V v v I 



" l = 4 ^ 4 /^(- r (^) 



i • def 

being p = 



. n 2 

X R (m/A) 



Here the parameter p can be computed by means of the Eq. (|27p . In the limit A — > oo it 

A«(l/A) 

approaches the unity regardless of the actual value for rh: for instance, we find p E [1.019, 1.072] for rh G [2, 12] (at 
A ~ 4 x 10 18 ). 

Implementing the Eq. (j2"9")l is straighforward and computationally not demanding. In spite of its simplicity, it clearly 
exhibits the iterative solution scheme, following the spirit of the typical condensed matter approach to the solution 
of the Dyson-Schwinger equations original sets (|1VI22[) . 

We find rh = 8.88 ± 0.10, where the error comes from the parameter p ranging in the aforementioned interval. 
Our finding is in agreement with the value of rh = ~ 8.88577 predicted by the classically scale invariant (CSI) 

theory and computed on the lattice [44j|. Indeed, in CSI theory an effective potential is computed starting from the 
cuclidcan action of a massless self-interacting scalar field and integrating out the field fluctuations around its vacuum 
expectation value (VEV). Different schemes are possible: in Ref. [Hj], for instance, two different renormalization 
constants are introduced for the field VEV and for the fluctuations. Another possibility is to add a mass term in the 
euclidean action and to describe the broken-to-symmetric phase as a phase transition. In order to account for different 
phcnomcnological approaches, the effective potential depends on a parameter ranging between 1 (CSI case) and 3 
(classical quartic potential), while the value 2 is a phase transition signature. The effective potential is then expanded 
up to the third power around its minimum and the three coefficients of the expansion, related to the abovementioned 
parameter, are fitted by the lattice simulation data. The lattice data fitted central value is 1.14, although the resulting 
uncertainty is large (see Ref. [44| and references therein for details). While the latter statement signals the need for 
more statistics in lattice calculations, the former result points at the CSI theory as the preferred scenario. 

In conclusion, we believe that the Dyson-Schwingcr equations in QED, together with the Green-Ward-Takahashi 
identity, we investigate in this paper, are equivalent to the analogous set of integral equations studied in condensed 
matter, often referred to as the Hedin's equations. They account for the self-consistency of the method, corresponding 
to the truncation of the perturbative/diagrammatic expansion in particle physics. Within the scheme proposed, the 
approach seems to be capable of producing the non-perturbative solution of the self-interacting scalar field theory, 
as suggested in the CSI theory, within theoretical errors. This method could be straighforwardly applied to other 
abelian (as well as non-abelian) theories, for instance the self-interacting fermion field theory, as an effective theory for 
describing the gap opening (mass generation) mechanism in condensed matter. We only mention here a result for the 
graphene, a system composed by sp2-hybridized carbon atoms placed on a plane. Such a system is interesting both on 
experimental and theoretical side, due to its peculiar electronic properties (see, for instance, Ref. [45| and references 
therein). While from the former point of view the gap opening is still under investigation, it is not theoretically 
understood yet whether or not its spatial dimensionality is two (as assumed in all the theoretical works) or three 
(as suggested by the role played by the off plane p z orbitals). Both issues (gap and dimensions) could be addressed 
by a Dyson-Schwingcr approach together with a similar mass operator ansatz: indeed the mass/gap equation, as it 
appears in the Eq. (|2"5"j) . actually depends on the space-time dimension (38|. Therefore, if a gap opens, it will depend 
on the actual spatial dimension (probably between two and three). The details of this work will be given elsewhere 
together with the results. 
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